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Abstract 

We derive an expression for the expansion of outgoing null geodesies in 
spherical dust collapse and compute the limiting value of the expansion 
in the approach to singularity formation. An analogous expression 
is derived for the spherical collapse of a general form of matter. We 
argue on the basis of these results that the covered as well as the naked 
singularity solutions arising in spherical dust collapse are stable under 
small changes in the equation of state. 



1 Introduction 



The spherical gravitational collapse of matter obeying an energy condition 
has been studied for various equations of state, including dust, perfect fluids 
and scalar fields. It is known that both black-hole and naked singularity 
solutions can result as the end-state of spherical collapse, though the issue 
of the stability of these solutions remains an open one. 

The expansion of a congruence of null geodesies, usually denoted by 9, 
can play a useful role in our understanding of the nature of the singularity 
that can form in gravitational collapse. With this in mind, we derive an 
expression for 9 for spherical dust collapse and also for spherical collapse of a 
general form of matter. Using these expressions we argue for the stability of 
both the covered and naked singularity solutions under small perturbations 
of the equation of state. 



2 Geodesic Expansion in the Dust Model 

In comoving coordinates the Tolman-Bondi metric, which describes spherical 
dust collapse, is given by 

ds 2 = dt 2 - e"dr 2 - R 2 dfl 2 . (1) 

The field equations are 
and 

R 2 = ^ + f- ( 3 ) 

The energy density p(t, r) is the only non-vanishing component of the energy- 
momentum tensor. The function /(r) above is defined by the relation 

(4) 
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and the function F(r) is twice the mass within the comoving coordinate r. 
Dot and prime denote differentiation with respect to the comoving coordi- 
nates t and r respectively. We take R to be negative as we are considering a 
collapsing cloud. 

Consider a congruence of outgoing radial null geodesies in this space-time, 
having the tangent vector (K f , K r , 0, 0), where K* = dt/dk and K r = dr /dk. 
The geodesic expansion 6 is given by 



which gives 



-g dx % 



dK l dK r K l , , x. K r f 



dt dr- 
ln order to compute the sum 



we proceed by noting that 

dK r 
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and similarly 
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dk dk dt dk dr 
Dividing the first of these two relations by dr/dk and the second by dt/dk, 
and after adding the two equations we get 



dK l dK T 

nr + — 



dr 



1 dK r 1 dK l 
+ 



K r dk K* dk 



4 4 



(10) 



In arriving at this equality we have used the fact that for outgoing radial null 
geodesies 



*L = * = e ^ 2 
K r dr 

Further, it follows from the geodesic equation that 



(11) 
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dk 
and 



dK r 

-V r rr {K r f - 2V r tr K t K r , (12) 



<UXl -Kri^f. (13) 



dk 

Together, these relations give the desired expression for 9 




(14) 



Let us now examine the relation that 9 bears with the formation or oth- 
erwise of a naked singularity in spherical dust collapse. As is well-known, the 
spherical collapse of a dust cloud results in the formation of a shell-focussing 
curvature singularity which lies on the curve R(t, r) = 0. (We assume that 
the initial density monotonically decreases away from the center, so that no 
shell-crossing singularities form). The central singularity, R(to,0) = 0, is 
known to be (at least locally) naked for some initial data and covered for 
other initial data, whereas the non-central singularity R(t, r ^ 0) = is 
necessarily covered. Consider first the simpler case of the non-central singu- 
larity; here the ratio F/R goes to infinity as the singularity is approached, 
while all other quantities in the expression for 9 are finite. As a result, 9 goes 
to negative infinity, which is consistent with the fact that the singularity is 
covered and not visible. 

The peculiar nature of the central singularity has been analysed in various 
previous papers [|IJ; in particular we draw attention to our analysis in [g]. 
There we check for the possible occurrence of a central naked singularity as 
follows. We define a quantity X = R/u, where u = r a , and a > 1 is defined 
such that R' /r a ~ l is a unique finite quantity in the limit of approach to r = 0. 
Then it is evident that in the approach to the singularity we can write the 
limiting value X of X as 



X = hm A = lim-=lim^ = lim-^ T (l-,/Z±M). (i 5 ) 

The quantities on the right of the final equality are written as functions of 
X and r, that is, X is used as a variable instead of t. The occurrence or 
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otherwise of a naked singularity depends on whether or not this equation has 
a positive real root. More explicitly, by calculating the limiting value of Bl 
it is shown that the singularity is naked if the equation 



admits a positive root for X. The quantity 7 is a known positive function of 
the initial data, and A(r) = F/r a . It is seen from the above equation that 
whenever the singularity is naked the expression inside the second bracket 
is positive definite. Further, it is known that whenever the singularity is 
covered, this same expression is negative in the limit. 

The connection with the limiting behaviour of 9 is seen as follows. By 
comparing the expression for 9 with Eqn. (|16|) above, we see that the limiting 
value of 9 is simply 2a /k if the singularity is naked. (We note that near 
r = 0, K r = dr/dk w r/k). Hence 9 goes to positive infinity at the naked 
singularity. On the other hand, whenever the singularity is covered, 9 goes to 
negative infinity. This kind of behaviour of 9 at a covered/naked singularity 
is of course to be expected. We merely intend to show that the occurrence of 
a central naked singularity is a result of the fact that for certain initial data 
the geodesic expansion for outgoing null geodesies continues to be positive 
even as the singularity is approached. 

We now attempt to ask how the information about the formation or 
otherwise of a naked central singularity is contained in the initial distribution 
for 9. What we do know of course is that for some initial density and velocity 
distributions the resulting singularity is naked. We can expect that somehow 
this must reflect in the initial profile for 9. As for the actual value of 9 
at any given point, that is naturally positive initially, for a congruence of 
outgoing geodesies - by itself this carries no information about the nature of 
the resulting singularity. Hence, we need to look at the initial distribution of 
9(r) and how it changes from one point to another. Consider the expression 
(|14D at a general time t for the marginally bound case (/ = 0). We recall 
that the solution for the area radius is 



where the scaling R — r at the initial epoch t = is assumed. Let us also 




(16) 




(17) 
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recall from that the series expansion for F(r) near r = is written as 



F(r) = F r s + F q r q+S + ... (18) 

where it is assumed that the first non- vanishing derivative in an expansion for 
the density near r = is the gth one. Using this, we can write the expression 
( |HD for 9 near r = 0, to leading order, as 



2Rf_ I F 3/ V 
R [ [l-l^t-^rn 



* 1/3 " 
o "" 1 



W,r) = — \l-\ - , ^ ^- 1 \K r . (19) 

4F ' 



We concentrate on the behaviour of the expression inside the second bracket, 

/F s3/2 F 3/2 3 

^(t,r)=(-) = g= ' . (20) 

This expression already contains information as to whether a naked sin- 
gularity will result or not, as we now elaborate. The term which carries 
information about the inhomogeneity in the distribution is the last term in 
the denominator: ( — jj^r q i). The remaining terms in i/)(t,r) are exactly 
as they would be for a homogeneous cloud. In the approach to the time of 
formation of the central singularity, which is equal to 2/3a/Fo, the inhomo- 
geneous term causes if)(t, r) to go to zero or infinity depending on whether q 
is less than three or greater than three. As a result, 9 goes to a positive limit 
if q is less than three, and to a negative limit if q is greater than three. This 
is consistent with the fact that the resulting singularity is naked for q < 3 
and covered for q > 3. The case q = 3 is however described satisfactorily 
by if; only to an extent: 9 is negative in the limit for ( = F^/F^ 2 < —4/3, 
whereas it is known that a naked singularity results only for ( < —25.9904. 

How can we look at the initial data for 9 and decide whether or not 
a naked singularity will result? We propose to look at the expression for 
l/i/)(t, r), which is 



if){t,r) \F) F 3/ V 1 ] 

If we look at the inhomogeneous contribution, given by the second term, at 
an epoch t just after the start of evolution, we notice that this contribution 
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diverges at r = 0, if q < 3 and goes to zero for q > 3. (We do not consider the 
transition case q — 3). In fact this continues to be so at all times t. Hence 
one could suggest that if the inhomogeneous contribution initially diverges 
at the center the singularity will be naked, and if it converges to zero, the 
singularity will be covered. Physically this means that at the initial epoch, 
the inhomogeneity can cause the geodesic expansion to either decrease or 
increase as one moves away from the center. The former case results in a 
naked singularity, and the latter in a covered singularity. The same set of 
arguments also hold for the non- marginally bound case (/ ^ 0). 

3 Geodesic expansion in the general spherical 

case 

We now show that the expression for 9 in the case of general spherical collapse 
can be cast in exactly the same form as for dust collapse. To begin with, 
we write the Einstein equations for the general spherical case in comoving 
coordinates. The metric, in comoving coordinates, is 



and the energy-momentum tensor is T\ = diag(p,p r ,pT,PT)- The field equa- 
tions for this system are 



ds 2 = e a dt 2 - e w dr 2 - R 2 dQ 2 



(22) 



F' 



(23) 



P = 



R 2 R r 



F — — 



p r R 2 R, 



(24) 



2p' r 4R' 



(PT-Pr), 



(25) 



P+Pr R(p + Pr) 



U = 



2p 4R(p + p T ) 



(26) 



p+p r R(p + p r ) ' 



and 




(27) 
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The function F(t, r) is equal to twice the mass inside the comoving coordinate 
r, and as in the dust case, the function f(t,r) is defined by the relation @. 
The difference from the dust case is that now the functions F and / depend 
on time as well. (It is only for the dust equation of state that both these 
functions are time-independent). When one considers the special case of 
dust, the above five equations behave as follows. Eqn. ( p3|) holds as such, 
while (^7|) reduces to @. Equations (|24j) and (|25| ) are trivially satisfied, 
whereas fl2~6|) reduces to an identity. 

The calculation of 9 proceeds exactly as in the case of dust, and we get 
again the expression in Eqn. (||). In order to calculate the sum (|7|) we use 
the relations © and @ and we note that now for outgoing null geodesies 
we have 

K l dt 



— = — = e ^~ u )/ 2 
K r dr 

Hence we get, by proceeding as in the dust case, 



(28) 



dK* 8K r 



dr 



1 dK r 1 dK* 
+ 



K r dk K* dk 



+ 



4 



4 



(29) 



By substituting for dK r jdk and dK l jdk from the geodesic equation we again 
find that 9 is given by 



6{t,r) 



2R' 



i f + F/R 
1 + / 



K r 



(30) 



which is the same expression as in the dust case, except that now F and / 
are functions of t as well as r. We have also utilised the field equation ([27]) 
which in fact is the only field equation used in writing the above expression 
for 6. 



4 Discussion 

We note that 9 for dust collapse has a peculiar behaviour in the limit of 
approach to the singularity curve. As regards the points other than the 
center (i.e. those with r ^ 0), they first get trapped (i.e. 9 becomes zero), 
and then singular. This is consistent with the singularity theorems - the 
formation of a trapped surface at a given r is followed by the formation of 
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a singularity for this value of r. However, r = is a very special point: 
the apparent horizon curve begins to form at r = simultaneously with 
the occurrence of the singularity. As a result, it would appear that one 
could not directly appeal to a singularity theorem to predict the formation 
of the central singularity, even though an explicit calculation shows that a 
singularity does form at r = 0. Similarly, one cannot predict a priori that 
the central singularity will necessarily be covered, and explicit calculation 
actually shows it to be otherwise. 

We would like to argue here, judging by the expression for 9, that both 
the black-hole and naked singularity solutions arising in spherical dust col- 
lapse are stable under sufficiently small changes of equation of state, for a 
fixed initial data. This is expected to hold for all initial data, except near 
the transition region ( i.e. the case q = 3) in the above discussion, and the 
reasoning is as follows. Consider an initial density and velocity distribution 
in spherical dust collapse which results in a covered singularity. As discussed 
above, this covered singularity results because for this initial data, the ratio 
Fj R goes to positive infinity in the limit of approach to the singularity. We 
have noticed, in Eqn. (|2lD , that the initial data already contains the infor- 
mation that F/R will behave in this particular manner in the approach to 
the singularity, thereby making it a covered singularity. Now, if we keep the 
initial density and velocity distribution fixed, and make a sufficiently small 
change in the equation of state, the initial behaviour of F j R will be arbitrar- 
ily close to that in the dust case. This is to be expected from the stability 
of solutions of Einstein equations under a small change in the equation of 
state. Hence the evolution will still be such that F/R will again go to pos- 
itive infinity as the singularity is approached, which makes 9 negative, and 
hence the singularity continues to be covered. (The function f(r) plays an 
insignificant role in this argument: one only has to make the plausible de- 
mand that /(0) be finite during the evolution). A similar argument ensures 
that the dust naked singularity solutions will be stable under perturbation 
of the equation of state (in this case F/R goes to zero in the limit and 9 is 
positive). It is only in the transition region (q = 3) that both the naked and 
covered solutions could be unstable: one kind of solution could go into the 
other kind. 

From an argument of this nature it also follows that if for any given 
equation of state some initial data lead to a covered (naked) singularity, the 
nature of the singularity will be left unchanged if the equation of state is 
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perturbed while keeping the initial conditions fixed. The same conclusion 
has also been arrived at earlier in || by independent means. It appears as 
if both the covered and naked solutions arising in spherical collapse must be 
treated on the same footing, in so far as stability under change of equation 
of state is concerned. 



References 

[1] D. M. Eardley and L. Smarr, Phys. Rev. D 19 (1979) 2239; D. 
Christodoulou, Commun. Math. Phys. 93 (1984) 171; R. P. A. C. New- 
man, Class. Quantum Grav. 3 (1986) 527; I. H. Dwivedi and P. S. Joshi, 
Class. Quantum Grav. 9 (1992) L69; P. S. Joshi and I. H. Dwivedi, Phys. 
Rev. D47 (1993) 5357. 

[2] T. P. Singh and P. S. Joshi, Class. Quantum Gravity 13 (1996) 559. 

[3] I. H. Dwivedi and P. S. Joshi, Commun. Math. Phys. 166 (1994) 117. 



10 



